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We investigate properties of ultralong-range polyatomic molecules formed with a Rb Rydberg 
atom and several ground-state atoms whose distance from the Rydberg atom is of the order of n 2 ao, 
where n is the principle quantum number of the Rydberg electron. In particular, we put emphasis 
on the splitting of the energy levels, and elucidate the nature of the splitting via the construction 
of symmetry-adapted orbitals. 



I. INTRODUCTION 

The recent advancement of ultracold physics has made 
possible the study of many interesting phenomena, rang- 
ing from the formation of molecules in a Bose-Einstein 
Condensate (BEC) 0, 0J1. to correlation effects in ul- 
tracold neutral plasmas [il H, where in the former case 
the atoms are cooled to temperature in the nano-Kelvin 
range. Combined with narrow bandwidth lasers and high 
resolution spectroscopy p , new phenomena can be stud- 
ied involving high- lying Rydberg states which have a nar- 
row spacing in energy of the order of 10 GHz. One 
such example is the theoretical prediction of the for- 
mation of ultralong-range dimers by a Rydberg atom 
and a nearby ground-state atom 0], the so-called "trilo- 
bite" molecules. The potential well supporting the vibra- 
tional bound states is extremely weak compared to typi- 
cal ground-state molecules. The depth of well is ~ 15GHz 
for n=30, where n is the principle quantum number of 
the Rydberg atom, and scales as n~ 3 . The long-range na- 
ture of such system, bound at the equilibrium distance of 
the order of 10 3 a.u., is rather unusual as well as the oscil- 
latory feature at the bottom of the potential well. This 
class of ultralong-range dimers should be distinguished 
from another kind, due to the Rydberg-Rydberg interac- 
tion @, whose molecular resonance was observed in the 
experiments |oL Hoj|. 

In this paper we address the question if more than 
one ground-state atom can form, together with the Ryd- 
berg atom, a polyatomic molecule. We find, that shape 
and symmetry of such polyatomic molecules follow a sys- 
tematics which is well understandable on the one hand 
side in terms of trilobite building blocks, that is, linear 
combinations of wavefunctions which describe the Ry- 
dberg atom and a single ground-state atom. From a 
more global perspective, these polyatomic molecules can 
be classified according to irreducible representations re- 
flecting their symmetry properties, as it is well known in 
quantum chemistry. 

So far, trilobites have not been identified experimen- 
tally. As the experimentally achieved density of ultracold 
atomic ensembles increases [ll|, so does the likelyhood 



of detecting such molecular species. However, there are 
still open questions concerning the realizability of such a 
molecule under current experimental conditions. For ex- 
ample, the Rydberg electron is very likely to interact with 
multiple ground-state perturbers, and it is not so clear 
at first sight what the role of the Rydberg states is. Such 
scenarios will also concern the proposed "dipole block- 
ade" scheme for quantum information processing 
The characterization of polyatomic molecules involving 
one Rydberg and several ground-state atoms helps to un- 
derstand the possible role of such ground-state perturber 
better. 

Technically speaking, we extend the calculations done 
by Greene, Dickinson and Sadeghpour 7] to include 
multiple ground-state atoms using the Fermi pseudo- 
potential treatment. More sophisticated methods exist 
[3 [li]], but the results do not differ much, whilst the 
qualitative features are certainly captured which suffice 
for the purpose of the present article. Using the Fermi 
pseudo-potential also allows one to calculate easily a large 
system. 

We investigate, in particular, the effect of placing the 
perturbers in a structured environment on the splitting of 
the adiabatic energy levels of the molecular system. We 
will use group theory to obtain the total wavefunction of 
the system in the framework of first-order perturbation 
theory via the construction of the symmetry- adapted or- 
bitals. For the sake of clarity we restrict the investigation 
to atoms all lying in a plane (whereas the Rydberg elec- 
tron of course lives in the physical 3D space). Atomic 
units are used unless stated otherwise. 



II. THE HAMILTONIAN 

First, we consider a ground-state atom with label i 
located at distance Ri from the Rydberg atom. The 
ground-state atom influences the electron by its polar- 
ization field, which has the form —a/2r , where a is the 
atomic polarizability. For Rb atoms, the experimentally 
determined value is a = 319.2 [l5|. To a good approxi- 
mation, the potential - extremely short-ranged with re- 
spect to the extension of a Rydberg electronic wavefunc- 
tion - can be mimiced by the Fermi pseudo-potential, 
namely [l6l |. 
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V l = 2^L[k l ]8{r-R l ). 



(1) 



2 



where L[ki] = — (tan5 s )/fcj is the s-wave energy- 
dependent scattering length of the Rydberg electron col- 
liding with a neutral atom; kf/2 = —l/2n 2 + I/B4 its 
kinetic energy, and r its distance from the mother ion. 
The s-wave phase shift S s can be calculated using the 
modified effective range theory by O'Malley et al 
and the zero-energy scattering length for triplet s-wave 
calculated by Bahrim et al [18j| , The singlet scattering 
length is much smaller, and hence we do not expect it to 
influence the phenomena discussed here. 

In the case of Rb atoms, the quantum defect is negli- 
gible for high /-states (I ^3), which are quasi-denergate. 
They are therefore well represented by hydrogenic wave- 
functions. The low-/ states, on the other hand, split away 
from the n-manifold and do not interact with the high- 
/ states provided that the energy-dependent scattering 
length is sufficiently small. The high-/ states are also 
more interesting, because the Hilbert space is larger in 
this case, which produces a more flexible system, i.e. the 
different eigenstates are allowed to interfere with each 
other. This results in a highly-polarizable complex. In 
this paper, we consider only the high-/ class. 

The total Hamiltonian of the Rydberg electron inter- 
acting with its mother ion and N ground-state atoms can 
be written as 



H = H + V N , 
where the ionic Hamiltonian is 

Ho = P -~ 1 -. 

2 r 



(2) 



(3) 



The potential Vn is the iV-fold sum over the interaction 
Eq. (fT]) with all ground-state atoms, i.e., 



N 



V N = 2n^2L[ki]6(r-Ri), 



(4) 



where i labels the i-th ground-state atom, and N is the 
total number of ground-state atoms. 

The effect of p-wave electron scattering was neglected 
in our calculation, however, in the dimer case it was stud- 
ied by Hamilton et al [19( and Khuskivadze et al [141 ] . 



III. DETERMINATION OF ADIABATIC 
ENERGY CURVES 
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FIG. 1: Adiabatic energy curves for the linear diatomic (dot- 
ted line) and triatomic molecule (solid lines) as a function of 
distance R between Rydberg and ground-state atoms, see also 
the sketch for the triatomic molecule with the Rydberg atom 
in the middle. The solid lines are the results from diagonal- 
ization, and the points (+) and (o) are the enegry expectation 
values calculated from the symmetry-adapted orbitals Eq. © 
and (0, respectively. The E = lies at the energy of the 
n — 30 manifold. 



using two methods: (i) the direct diagonalisation of H; 
and (ii) the projection operator method [20| to construct 
symmetry-adapted orbitals and determine the BO curves 
from standard perturbation theory. Both methods are 
accurate, but the latter gives a deeper insight into the 
quantum mechanical properties such as the energy de- 
generacy. 

The Fermi pseudo-potential is usually valid for n >25- 
30. For smaller principle quantum numbers, the scat- 
tering of e~ + Rb and the polarization of the neutral 
perturber by the Rydberg atomic core are not indepen- 
dent of each other [2l|. Here, we present calculations for 
n = 30, which also allows us to compare our results di- 
rectly with that previously obtained for the Rb2 dimer 

0.Q3. 

The eigenvalues from the unperturbed Hamiltonian, 
Ho, yields simply the hydrogenic energy Eq = —l/2n 2 , 
so for convenience we set this to be zero throughout this 
article. 



Under the adiabatic approximation, the energy levels 
calculated from the Hamiltonian H (including the per- 
turbation) give, automatically, the electronic structure of 
the molecular system involving N + 1 atoms. We deter- 
mine a cut through the Born Oppenheimer (BO) poten- 
tial surface for systems with TV = 2, 3, and 4 ground state 
atoms which are uniformly placed on a circle with radius 
R centered at the Rydberg core. These three cases cor- 
respond to linear, triangular, and square geometries, re- 
spectively. The cut we choose corresponds to the breath- 
ing mode, i.e. R is varied. We calculate the BO curves 



IV. COLLINEAR TRIATOMIC MOLECULE 
(N=2) 

Consider two ground-state atoms (A and B) placed on 
either side of a Rydberg atom with distances R forming 
a collinear triatomic molecule. This configuration cor- 
responds to N—2, and the numerical result of the BO 
curves are plotted in Fig. [1] 

In order to understand the splitting of the energy lev- 
els, we use the perturbed state \ip n (r)) when only one 
of the two ground-state atoms is present as the building 
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block for constructing the total electronic wavefunction. 
The perturbed state can be can be explicitly written as 



iv>»(o>=E#u$)iM*)>. 



(5) 



where the index q runs over all the degenerate states 
which includes all Vs and m's with I > 3. We call this 
wavefunction the "trilobite" wavefunction since it pro- 
duces the probability density like that drawn in Ref. 
The two wavefunctions, which clearly satisfy the parity 
of the collincar geometry, arc 



and 



|V+(f)) = |^(r)) + |^(r)>, 



\^(r)) = \^(r))-\^(r)), 



(6) 



(7) 



where the superscript A and B are the labels of the 
ground-state atoms; |V'n(^*)} an d IV'n(^)) are t ne trilo- 
bite wavefunctions when only atom A or B is present. 
By choosing the projection axis z so that it aligns with 
the internuclear axis, the only degenerate states that con- 
tribute are those with non-zero value along z. They are 
in this case the states with m = 0. 

Using the above ansatz to calculate the expectation 
value (Vjv) yields immediately two energies E+ and 
which are distinguished by their parities: 



n-1 



E ± (R) = L[k}J2(2l + l)\uni(R)\ 1 



1=3 



I = even for E + , 
I = odd for E_, 



(8) 



where L[k] = L[kA\ — £[fcs], R — \Ra \ — \Rb \ and u n i is 
the radial part of the hydrogenic wavefunction <fi n i(R) = 
U n i(R)Yi m (6,ip). 

Hence, we see that, with the inclusion of the second 
perturber, two curves split away from the n-manifold, 
with one corresponding to gerade and the other to 
ungerade symmetry. They both converge at large dis- 
tance to the curve when only one ground-state atom is 
present. However, they split from each other approxi- 
mately within R/r n < 1 (with r n — n 2 — 900), which 
can be seen from the sum of the probability densities of 
even-Z and odd-l states. They differ when the overlap 
(ipn\ipn) i s not exponentially small, see Fig. [2] This fea- 
ture is general for all principle quantum numbers. The 
additional splitting also suggests that the system can be 
more stable with the inclusion of more than one neutral 
perturber, a situation we will investigate in more detail 
in Section Ivl 

In Fig. [31 we show the contour plot of the probabil- 
ity density of the diatomic and triatomic systems at the 
interatomic distance R m — 879, which corresponds to 
the deepest potential energy. The special minimum con- 
figuration can be most clearly identified by means of the 
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FIG. 2: Sum of probability densities 2TJ ; |0 n j(J?)| for states 
with even and odd angular momentum components / as a 
function of the radial distance. 




FIG. 3: Cut of the electronic probability density along the 
internuclear axis at the deepest point of the potential well, 
R = 879. The black contour lines show the probibility density 
of the trilobite (diatomic) wavefunction, and the background 
coloured plot is for the N = 2 (collinear triatomic) configura- 
tion. The ground-state atoms are located at [x, z) ±879) 
and the Rydberg atom at (0,0). The two white solid lines 
show the classical Kepler ellipses. 



classical Kepler orbits along which the two trilobite states 
of the molecule are scarred (see also Ref. [23]). Each Ke- 
pler ellipse has one ground-state atom in one of its foci 
and touches the other ground-state atom. 
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FIG. 4: Illustration of the geometry of the (a) C3,„ and (b) 
Ga,v configurations. The dashed lines are the planes of reflec- 
tion. The shaded and the solid circles are the ground-state 
and the Rydberg atoms respectively. 



V. PLANAR POLYATOMIC MOLECULES OF 
TRIANGULAR AND QUADRATIC SHAPE 

(N=3,4) 

For N > 3, we choose z to be perpendicular to the 
plane containing the atoms. Figure 0] illustrates the spa- 
cial geometry of the complexes. The degenerate states 
that contribute, i.e., the index q in Eq. (J5]), are the states 
with I + m being an even integer, exclusive of I = 0, 1, 2 
states. The results are shown in Fig. [5] 

For the N = 3 case, three energy eigenvalues split away 
from the n-manifold, and as in the case of N = 2, be- 
yond r n they converge to the BO curve for the dimcr 
case. Note, however, that two of the energy levels are 
degenerate at any given distance within r„, indicating 
that there are additional symmetries preserved under the 
perturbation of the ground-state atoms. To elucidate 
these symmetries, we again use the trilobite state as ba- 
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FIG. 5: The three adiabatic energy curves Qi for a triangular 
configuration as a function of Rydberg-ground state atomic 
distance R (see Fig. [4^,) . The coding of the data is analogous 
to Fig. [TJ 



TABLE I: Character tables of the irreducible representations 
of (a) Cz,v and (b) Ci lV [25|. The labeling follows the conven- 
tional rules. 



(a) 



(b) 
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sis functions to construct the relevant symmetry-adapted 
orbitals. But unlike in the previous collinear configura- 
tion, where taking account of the parity as the relevant 
symmetry is intuitive, we have to use a systematic ap- 
proach for N=3 or larger. A method which has been 
used extensively to find the symmetry-adapted orbitals is 
the projection operator method, which is expressed math- 
ematically as [20] 



(9) 



where Or is the operator for a particular symmetry op- 
eration i?, e.g., rotation or reflection, etc.. The lj and 
Xj are, respectively, the dimension and the character of 
the j-th irreducible representation (irrep) of the symme- 
try group, to which the system belongs, while h is the 
order of the group. The sum extends over all symmetry 
operations in the group. This equation allows us to find 
the symmetry-adapted functions {<^} from any original 
basis set {/j}. The general proof of this proceedure can 
be found, for example, in section 6.6 of Ref. (2~i ]. 

The name of this procedure originates from the fact 
that the pre-factor in front of / in the above equation 
can be viewed as a projection operator that projects the 
basis set {/$} into a basis set {<?i} that diagonalizes the 
Hamiltonian matrix. In other words, g and / in Eq. §§§ 
are vectors, and the operator 



Pr = 



R 



(10) 



can be represented by a unitary matrix. 

The configurations of 7V=3 and 4 correspond to the 
symmetry groups C^.y and C^ Vl and their character ta- 
bles of the irreducible representations are shown in Ta- 
ble D 

We find that in the N=3 case, the representation of the 
symmetry operations using the trilobite state as the basis 
set contains only two of the total three irreducible repre- 
sentations, Ti and (see Table 1(a) and Appendix IA II 



for details), which are one- and two-dimensional respec- 
tively. The symmetry-adapted orbitals constructed will 
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FIG. 6: Same as Fig. [5] but for a square geometry with a 
ground-state atom at each corner (see Fig. 2b) with four en- 
ergy curves gi, see text. 



then, according to the fundamental theory of quantum 
mechanics, consist of a non-degenerate and two degen- 
erate states. They are, repectively, g%, gi and 5(3 shown 
explicitly in Eq. (|A"5)) , (|A6a[) and (|A"7|) in Appendix I A 11 
The energy expectation values, 



(V N (R)) gi = (9i\V N (R)\ 9i }, 



(11) 



calculated using the ^-functions i = 1,2,3 are plotted 
in Fig. 03 As expected, the two curves belonging to T 3 
overlap with each other at all distances R. Note that 
applying Pp 2 onto any of the basis functions produces 
zero, which is a general feature when the irrep is not 
contained in the overall representation. 

The same analysis for the C^ v symmetry reveals that 
the overall representation contains Fi, T3, and T$ (see 
Table. 1(b)), with the first being one-dimensional, and 
the second and the third two-dimensional (see Appendix 
IA 21 for detail). Hence, there are two sets of doubly- 
degenerate BO curves and a non-degenerate one. Ap- 
plying Eq. ^ , we obtain the symmetry-adapted orbitals 
gi, i = 1, 2, 3, 4 as shown in Eq. (|A11|) - (|A14|) in Appendix 
IA 21 The adiabatic energy levels from the analytical and 
numerical results are plotted in Fig. [S] Again, the graph 
shows a perfect agreement between the two results. 



VI. CONCLUSION 

We have used the Fermi pseudo-potential to model the 
effect neutral Rb ground-state atoms have on a Ryd- 
berg electron. Cuts through the resulting potential sur- 
face, adiabatic in the distance of the ground-state per- 
turbers from the ionic core of the Rydberg electron, have 
been calculated for different arrangements of the ground- 
state atoms which form planar polyatomic molecules. 
We found that more ground-state atoms lead to more 
deeply bound molecules compared to the original diatom 



as studied by Greene and coworkers Q. A systematic 
understanding of the structure and symmetry of such 
molecules can be gained by taking the trilobite (diatomic) 
wavefunction as a basic unit and constructing symmetry- 
adapted orbitals as demonstrated in Section IIVI and [V] 
For two ground-state atoms the procedure is relative sim- 
ple and intuitive, while three or more ground-state atoms 
require a systematic approach, such as the projection op- 
erator method. 

In the case where N is larger than the number of de- 
generate states <7max, our method of constructing the per- 
turbed wavefunction should still work, and will yield only 
gmax linearly independent states. 

The effect of p-wave electron scattering plays an impor- 
tant role especially in hydrogen- like atoms. Our calcula- 
tions can be extended into the case of higher partial wave 
scattering by using the appropriate pseudo-potentials for- 
mulated by Omont (22J. Previous studies [TJ, [l|| for 
dimers show that the potential curve of the p-wave scat- 
tering crosses the potential well at R « 1200, and could 
potentially destabilize the trilobite while also providing 
additional potential wells. It will be also interesting to 
see how the spatial arrangement of the atoms affect the 
energy of the system in this process, since the potential 
is now dependent on the gradient of the electronic wave- 
function in the 3D space. 

The present work is a first exploration of the possibil- 
ity to form polyatomic molecules from a single Rydberg 
atom and a number of ground-state atoms. We have only 
determined a cut (at equal distances of the ground-state 
atoms to the Rydberg atom) through the multidimen- 
sional potential surface which resembles the potential for 
the breathing mode of the molecule. Future analysis and 
realistic assessment of the quantitative features of such 
species must include the vibrational motion of the atoms. 



APPENDIX A: DERIVATION OF THE 
SYMMETRY-ADAPTED ORBITALS 

1. Planar Polyatomic Molecule with N=3 

The molecule formed in this configuration has the sym- 
metry of the point group Cq jV . Using trilobite states as 
the basis set to contruct the corresponding representa- 
tion, one obtains the following matrices, 



E = 



Or 




Or.x = 



O a ..„ = 




(Ala) 



(Alb) 



(Ale) 



6 



In the above notations, E is the identity; C\ denotes 
the rotation about z-axis by angle 2np/k] and the cr's 
are the relfections through the planes perpendicular to 



the plane of the atoms, as indicated in Fig. 4(a) 

It is clear that when one of the above operators, say 
O c i , is applied on the original vector, the result is a 120° 
rotation about z in the counter-clockwise direction, as- 
suming that z is pointing perpendicularly into the paper, 
i.e. 



(A2) 



From Eq. (|A1[) . the character of the representation in 
such basis set can be determined by taking the trace of 
each matrix, and they are summarized in the table below: 






1 














:) 




















E 2C 3 3a v 


Ttotal 


3 1 



Here, we have used Ttotai to denote the representation 
formed by the trilobite states. By inspecting the charac- 
ter table of the irreps of C 3jl , (Table 1(a)), one immedi- 



ately sees that the current representation is a direct sum 
of Ti and r 3 , namely, 

r to tai = ri©r 3 . (A3) 

Now, we determine the projection operators in each irrep 
by using Eq. (flT))) . The order of the group is h = 6, and 
the dimensions for Ti and r 3 are l\ = 1 and Z 3 = 2, 
respectively. Equation (|TUJ) then yields 

1 r. 
6 



Pi 



o E + o ( 



O 2 + o 



Pr 



Pr 



l r 

6 

l r 

6 



Oe + O nl + O c 



o 



o 



20 f 



Od + Or 



(A4a) 

Oa v3 ] , 

(A4b) 
(A4c) 



Since r to tai contains only Ti and T 3 , we need to apply 
only Eq. (|A4a|) and (|A4c[) to our basis set in order to 
obtain the symmetry-adapted orbitals. Acting the triv- 
ial operator Pp 1 on the trilobite wavefunction \4>n)i we 
obtain the first symmetry-adapted orbital 



P r J A =P Tl \^n) 



(A5) 

The same equations are obtained if one acts Pr 1 on \ip^) 
or \ipn) which are obviously linearly-dependent. How- 
ever, acting Pr 3 on > IV'n) an d IV'n) gives, repec- 
tively, 



|^}-|V£>] =91, 



PvM) = \Wi) 

PvM) = \Wl)-\€)-\<)\ 

P T Mn) = \[Z\€)-Wn)-Wn)] 



(A6a) 
(A6b) 
(A6c) 



Since T 3 is a three-dimensional irrep, two of the above 
equations can be combined, by subtracting Eq. (|A6b|) by 
(IA6c|) . giving 



53, 



(A7) 



so that finally we have three linearly-independent wave 
functions, which we call g\, gi and <? 3 . 



2. Planar Polyatomic Molecule with N=4 

Following the same procedure as in the case of N = 3, 
one finds the matrices of the symmetry operations in the 
point group Cn iV in the present basis set as, 



O f = 



O n i = 



o* vl 



o 



/iooo\ 

10 
10 
\0 1/ 

fo 1 o o\ 

10 
1 
\1 0/ 

(l o o o\ 

1 
10 
\0 1 0/ 

fo 1 o o\ 

10 
1 
\0 1 0/ 



O n i = 



O n s = 



Oa v 



O a 



fo i o\ 

1 
10 
\0 1 0/ 

fo l\ 
10 
10 

\0 1 oj 

fo 1 o\ 

10 
10 
\0 lj 

fo l\ 
10 
10 

\1 0) 



(A8a) 



(A8b) 



(A8c) 



(A8d) 



where the notations are as before, and the planes of re- 
flections are indicated in Fig. |4(b)| 

The character of r to tai can again be determined by 
taking the trace of each matrix above, and they are: 



E Ci 2C4 2(7^ la d 



total 



4 







Ag ain, from the character table of the irrep (Table 
|I(b)[ ), one finds that the representation r to tai is a direct 
sum of 



r to tai = ri©r 3 ©r 5 . 



(A9) 



Hence, we know that in this representation there are two 
one-dimensional and one two-dimensional irreps. Their 
corresponding projection operators can be obtained by 
applying Eq. (|10p . where in this case, h = 8, and I3 and 
I5 are 1, 1 and 2, repectively. Therefore, the projection 



7 



operators are 



following four linearly-independent equations: 



Pr, = l[OE + 6 ci +6 cl +6 c 



o n 



o„ 



1 



Pr, = -\0 E + O c i-0 



Or 



o n 



1 5 A 



O e - O c i 



(AlOa) 



(AlOb) 



(AlOc) 



The symmetry-adapted orbitals can then be obtained 
in a similar way as in Appendix IA 11 which yields the 



9i = 7 [lt A M) + lt B «) + \€(r)) + \^(f))} 



(All) 



93 = \ [Hn(r)) ~ 

.94 = \ [|Vf(r)) - \^{r))} 



(A12) 
(A13) 

(A14) 
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